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Abstract 

We develop a theory of Lagrangian reduction on loop groups for completely integrable 
systems after having exchanged the role of the space and time variables in the multi¬ 
time interpretation of integrable hierarchies. We then insert the Sobolev norm in the 
Lagrangian and derive a deformation of the corresponding hierarchies. The integrability 
of the deformed equations is altered and a notion of weak integrability is introduced. We 
implement this scheme in the AKNS and SO(3) hierarchies and obtain known and new 
equations. Among them we found two important equations, the Camassa-Holm equation, 
viewed as a deformation of the KdV equation, and a deformation of the NLS equation. 


1 Introduction 


The classification of integrable systems through hierarchies of commuting flow s such as th e 


AKNS hiermchy i s a well established theory which star t ed with Ablowitz et al. IQT.d 1974 ] 


Date et al. 1983l |: Flaschka et al. 1983a ]: Newell 1985 ]: Adler et al.l 2000 ] and encompasses 


almost all other notions of integrability, such as multi-Hamiltonian structures, Lax pairs, zero 
curvature relations (ZCR), r-functions, bi-linear equations and Pa inleve hierarchies. We refer 
the interested reader to the very complete book Scott et al. 2006| | and the references therein 
for more details on various other subjects in the theory of integrable systems. Recently, the 
discovery and the study of equations involving non-local dispersio n such as the Camass a -Holm 
equation op e ned a new area in i ntegrable systems. We ref e r to Camassa and Holm _ H)93l | ; 


Fokad 1995 ]: Fuchssteiner 1996 ]: Olver and Rosenau 1996 ]: Qiao 2007l | : 


Novikov 1200' 


iiffor 


some well studied equations of this type. Some of these equations are even physically rele¬ 
vant as higher approxima t ions of shallow water equation s. We refer to ICamassa and Holm 
1993]; Dullin et al. 2004 ]: Constantin and LannesI 2009l | for physical derivations of the CH 


equation. From this physical viewpoint, they are deformations of classical integrable systems 
or higher order approximations of more complete physical models. Despite these facts, it is 
well-known that the integrability of the deformed equation s is slightly different from t h eir clas¬ 


sical c ounterpart. They have non-local conservation laws (|Camassa and Holml 199311: 


2005(1), ZCRs without associated Zakharov-Shabat spectral problems ((Hone and Wand [2003 j 


Lenellsl 
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Constantin et al. 2006| | ). Perhaps these difficulties explain why a classification of these equa- 


tions, based on hierarchies such as the AKNS hierarchy, is still missing. Notice that recently 
Novikov 2009| | made a classihcation using an ansatz for the form of the equations and a test 


for their integrability, developed in Mikhailov and Novikov 2002l |. We will not follow their 


approach here because our aim is to understand each equation as a member o f a hie rarchy 


only defined with a Lie algebra in the sense of lAblowitz et ahl 1974i | ; iNewelll 1985l | . The 


key element that we will be using to develop such a theory is the fact that these deformed 
equations correspond to classical equations when the parameter a of the Helmholtz operator , 
or norm, is set to 0. For example, the Camassa-Holm equation Camassa and Holm 1993l | 
corresponds to KdV and the modified Camassa-Holm equation Fokas 1995 1: Qiaol ~ 2nr)7l | to 
mKdV. We will thus deform classical integrable hierarchies such that the deformed equations 
will be recovered and shown to correspond to a particular member of the original hierarchy. 

Following Date et al. 1983]; Flaschka et al. 1983a j: Newell 19851 ] we will use the loop 
group and multi-time interpretation of integrable hierarchies. The concept of multi-times is 
fundamental in this formul ation and makes sense of reduction procedures on the cotangent 
bundle of loop groups; see Pressley and Seeal 19861 ] for a detailed account on loop groups. 
In order to allow an equivalent Lagrangian formulation, we will extend these ideas by simply 
having a different interpretation of the multi-times. In the standard theory, the space variable 
is fixed and the flows of the hierarchy, or higher order symmetries, are spanned by the time 
variable. In the new Lagrangian interpretation, the time is fixed and the hierarchy is spanned 
by the space variable. Notice that our Lagrang ian theory is different from the pluri-Lagra ngian 
system s initiated in Lobb and Niihoff 2009 ] and further developed, for example, by SurisI 


20131] . 


In a second part, the usual 1? norm in the Lagrangian will be replaced by the norm 
and the corresponding deformation of the hierarchy computed. This use of the Sobolev norm 
is common to derive the CH equation as a geodesic m o tion on the group of diffeom or phisms 


of firs t Sobolev class; see for instance Misiolek 1998| ]: Holm and Marsden 20051 ] 


or 


Gnha 


20071 ] for yet another use. This procedure allows us to deform the entire classical hierarchies 
such as the AKNS hierarchy in order to recover the CH equation among others. The classical 


0 


Standard equation 

Deformed equation 

Limit —)■ oc 

sl(2) 

sl(2) 

(1,2) 

(1,3) 

NLS 112.261) 

KdV (12.231)/ mKdV (12.241) 
CmKdV (12.251) 

CH-WLS" (13.181) 

CH (13.131)/ mCH (13.151) 
CmCH (13.161) 

HS-NLS" (13.311) 
HS (13.301) 
mHS* (13.321) 

so(3) 

(1,2) 

(12.351)’' 

(13.251)’' 

- 

so(3) 

(1,3) 

(I2..39I)’' 

(13.281)’' 

(13:321)’' 


Table 1.1: Summary of the equations derived in this work using the hierarchy classification. 
The third column corresponds to classical equations such as NLS or KdV, the next column 
their deformations and the last one exposes a few limiting cases with o? ^ oo. We only 
considered the first two flows for the Lie algebra s[(2) and so(3) but other flow and Lie 
algebra could be derived in the same way. Asterisks indicate possibly new equations. 


integrable equations and their deformations can thus be classified though the Lie algebra g 
and the choice of space and time variables, or 2-dimensional slices Wj indexed by (i, j). This 
is summarized in the table 11.11 where the asterisks denote possible new equations. Almost 
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all the deformed equations are already known except the deformation of the NLS equation, 
which reads 


imt + Uxx + 2am{\u\‘^ — a'^\ux\‘^) = 0, m = u — a^Ux 


a = ±1. 


( 1 . 1 ) 


Indeed the weak integrability presented here does not guar antee its complet e integrability. De¬ 
spite the possible non-integrability, it has been shown in Arnandon 201(11 ] that the equation 
contains solitary waves and even peaked standing waves with almost elastic collisions. 
We want to mention that the CH - NLS equation is diff e rent f rom the generalized NLS equa- 
tion, first derived in FokasI 1995 ]: Olver and Rosenau 19961 ] and more recently studied by 


Lenells and FokasI 20091 ]. For the derivation of this equation they used the bi-Hamiltonian 


property of the NLS and CH equations to find an integrable extension of the NLS equa¬ 
tion, without using the Helmholtz operator in an intrinsic way. Other similar attempts for 
improving the NLS equ a tion, but without a sking the integrability question, were made by 
Colin and Lannes 20091 ]: Dumas et al. 20161 ] with an improved dispersion, also involving the 
Helmholtz operator. 


Structure of this work We will develop in section ITT] the Lagrangian d escription of inte- _ 

grable hierarchies with central extensions, similarly to the i2-matrix theory of ISemenov-Tvan-Shansku 


19831 ]. Th e Lagrangian reductio n theorem with central extension will be stated with the for¬ 


malism of iMarsden et al.l 20071 ]. It is worth mentioning that it seems to be the first time 


that this theorem is stated in this form, with this type of central extension. This new inter¬ 
pretation of the multi-times is then described in section 12.21 as well as how the corresponding 
Euler-Poincare or equivalent Lie-Poisson equations arise on two dimensional slices of the multi¬ 
times. On these slices, the time coordinate will be the usual dynamical coordinate and the 
space coordinate will be seen as the parameter of an infinite dimensional group or Lie algebra. 
The dynamics along the space coordinate will then be made non-trivial with the help of the 
derivative cocycle. Examples such as the AKNS hierarchy with SL{2) and another hierarchy 
with SO{3) will be shown in section (2^ After having set up the Lagrangian reduction, the 
deformation using the Sobolev norm is straightforward to implement in the Euler-Poincare 
equation. In section [3l the deformation of the hierarchy is derived and its integrability is 
investigated. As opposed to the classical case, where there is an equivalence between the 
Euler-Poincare equation and the associated isospectral problem, the Euler-Poincare equation 
cannot be directly interpreted as a ZCR. With the Eourier decomposition of the loop algebra 
elements, parts of the Euler-Poincare equation are trivially satisfied for the highest powers of 
the loop, or spectral parameter A, but are not valid in the deformed Euler-Poincare equation. 
We must therefore define a projection which removes these terms and makes sense of the 
projected Euler-Poincare equation, or projected ZCR. The corresponding PDEs will then be 
said to be weakly integrable if they satisfy a projected ZCR. 

This method allows us to deform all members of the AKNS hierarchy in order to recover equa¬ 
tions such as the dispersive Camassa- Holm equat ion Camassa and Holm 19931 ]: Dullin et al. 


2004 ] and the new CH-NLS equation Arnandon 2016| ]. This will be done in section 13.21 for 


AKNS hierarchy and then for the SO{3) hierarchy. 
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2 Lagrangian interpretation of integrable hierarchies 

2.1 Reductions with a central extension 


In this work we will consider a particular type of reduction by symmetry where the config¬ 
uration manifold is the group of symmetry itself. The corresponding reduction is called the 
Euler-Poincare or the Lie-Poisson r e ducti on, for respectively the Lagrangian or Hamiltonian 
mechanics; see iMarsden and Ratiul 1999l | for a complete treatment. The Lie group in this 


section will be infinite dimensional and of the form Map(M, G) where G is a Lie group. The 
M variable will be the space variable x of the 1+1 nonlinear PDEs that will be derived. The 
dynamics with respect to x will be made non-trivial by using a central extension with a 
derivative cocy cle. This system is different from usual 1+ 1 PDE s coming from a reduction by 
symmetry; see Ellis et al. 2010| |; Gav-Balmaz and Ratiu 2009| | for example. Indeed, in the 
standard theory, the dynamics on the space variable comes from an affine action of the group 
of symmetry on the advected quantities. Here, the dynamics arises from a central extension 
with a cocycle. 


2.1.1 Central extension 


We refer to IMarsden et al.l 2007l | for a complete treatment of group extensions in mechanics 
and we will only recall useful facts without proofs. A central extension := G x P of a group 
G by a vector space V is characterized by the action of G^ onto itself with a cocycle term in 
the extension of the group. A group two-cocycle is a bilinear map B{g,h) : G x G — ?• M which 
satisfies a cocycle identity such that the group action [g, a) ■ {h, b) = {gh, o + 6 + B{g, h)) is 
associative. The Lie algebra of G'^ is centrally extended by the tangent space of the vector 
space V. We will always use P = M and thus g'^ := g x M. The group cocycle drops to the 
Lie algebra by differentiation to give a Lie algebra cocycle c(^, ry) : g x g —)■ M which satisfies 
a cocycle identity such that the corresponding Lie bracket satisfies the Jacobi identity. The 
adjoint and coadjoint actions are given by 


ad(^,a)(ri,b) = l((,a),(?],b)] = and (2.1) 

ad(^,a)(+>™') = (ad|Ai + "rc(^,-),0), (2.2) 

where (^,a), (g,b) € and {g,m) G (g'^)*. We will also need the formulas for the inverse of 
a group element and for the tangent of the left translation 

= {g~^,-B{g~^,g)), {g,0)~^{g,0) = {g~^g,-D 2 B{g-^ ,g)), (2.3) 


where D 2 B stands for the derivative in the second slot of B. 

For the present theory there will be a space variable x. The dynamics along this variable is 
assumed to be smooth and will be given by the derivative cocycle B(g,h) = f gdxhdx. The 
corresponding Lie algebra cocycle is 

c{C,v) = j {C,dxg)dx, (2.4) 

where (•, •) is the Killing form on the semi-simple Lie algebra g. We will always consider 
semi-simple Lie algebras and periodic or vanishing boundary conditions. The main point is, 
as always, to identify g with g* and to be able to freely perform integrations by parts. 
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2.1.2 Lie-Poisson equations with a central extension 


The reduction procedure on the cotangent bundle of a Lie group leads to a Lie-Poisson equa¬ 
tion on the dual of the Lie algebra of this group. When using a central extension of the 
Lie group, the variable in the centre of the Lie algebra will always be a constant and thus a 
standard kinetic term can be taken in th e Hamiltonian. We refer to iMarsden et al.l 20071 ]: 
Garci'a-Naranio and Vankerschaverl 201.111 fo r more details of this construction. The theorem 
can now be stated; see iMarsden et al.l [2007t | for the proof. 


Theorem 1. Let be the central extension of the Lie algebra 0 with cocycle c : g x g 
and h : (g'^)* —>■ M 6e Hamiltonian function. The Lie-Poisson bracket is 


{F,G}{iL,a)) = (L, 


5F^ 5a 

Jl'Jl 


5F 5G\ 


and the Lie-Poisson equation is 


dtL = ad^L -F ac ( 

5L \ OL 


(2.5) 


( 2 . 6 ) 


These equations simplify, by using the derivative cocycle, the Killing form and a = 1, to 


{F,G}(L)= I 


L. 


Jl'Jl 


^ , ,5F ^ 5G\ ^ 

I \ 


e 7 

dtL — dx—= ad^L. (2.7) 

dL SL 


Notice that the form of the Lie-Poisson equation is the same as the usual zero curvature 
relation of integrable systems and is also the Lie-Poisson equa t ion used in the i7-matrix 
deriva tion of integrable system; see Semenov-Tvan-Shanskii 19831 ]: Blaszak and Szablikowski 
20091] . 


2.1.3 Euler-Poincare equations with central extension 

Provided that the Legendre transformation exists, the derivation of the corresponding Euler- 
Poincare equation is straightforward. However, in the integrable systems context, there is no 
Legendre transformation and the Euler-Poincare equation must directly be derived from the 
variational principle. We must therefore state the following Euler-Poincare reduction theorem. 

Theorem 2. Using the above definitions, the following statements are equivalent: 

(1) Hamilton’s variational principle with Lagrangian .if = J" .ifdx : G —)• M holds on G’^ 


5 / J^{g{t, x), g{t, x))dt = 6 // ^{g{t,x),g{t,x))dxdt = d 


for arbitrary variations 5g vanishing at the endpoints; 


(2) g{t) satisfies the Euler-Lagrange equations on G; 
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(3) the constrained variational principle 


5 


l{M{t, x))dt = S 


£{M{t,x))dxdt = 0 


holds on Q^, using variations of the form 


5(M,0) = (77+[M,r/],c(7?,M)), 


for arbitrary rj vanishing at the endpoints; 

(4) the Euler-Poincare equation with central extension holds, that is 


d 5i 
dt 6M 


1* \ 


( 2 . 8 ) 


(2.9) 


As for the Lie-Poisson equation, the derivative cocycle and the Killing form help simpli¬ 
fying the variation and the Euler-Poincare equation 

( /* \ X /? K /) 

dtr]+[M,ri\, j {r],dxM)dx\, 'diJM ^ ^ 

Only the proof with derivative cocycle and semisimple Lie group will be given as the general 
case is not of interest for this work. 


Proof. The equivalence of (1) and (2) comes from general theory of Hamilton’s principle. 
For (3), the reduced variations are computed using the action on the central extension of G 
provided by the cocycle B{g,h) = f gdxhdx. With the left trivialization formula (|2.3p . a left 
trivialized generic element reads 

(M,0) = to.0)-'te.0) = (s-^,-/s-‘4#<<.), 

and its variation decomposes as 

(5(M,0) = 5{g~^{g,Q)) = (^{g~^g),- j 5{g~^dxg)da^ , 


where the first slot gives the usual variation, namely 5M = ?)+ [M, rj\ for arbitrary g = g ^6g. 
The second term needs a little computation with integration by parts 

- j S{g~^dxg)dx = j [-g~^5gdx{g~^g) + g~^Sgg~^dxgg~^g - g~^dxgg~^Sg] dx 

= j [gdxM + gg~^dxgM - g~^dxgr] - g~^dxgMg] dx. 

Then, by noticing that 

0 = ^ y dx{g~^gg)dx = ^ Jdxg~^gg + g~^dxgr] + dxgdx = jg~^dxgr] + dxV, 
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the second term in the previous calculation vanishes with proper boundary conditions and 
similarly for the last two terms by using 

0 = y dxir]g~^gM)dx = 2r]g~^d^gM + d^{r]M) 

and ® ^ y ^^i 9 ~^ 9 Mr])dx = 2g~^dxgMr] + dx{Mg). 

We can now compute the Euler-Poincare equation with central extension from the variational 
principle and prove (4) 

5 y l{M)dt = yy + [^,»?],c(r?,M))^ dtdxdX 

5£ f 

— {ii + [M,r]])dtdx + J c{r],M)dt 

6£ \ 

{-dt + &(!%[)— + da:M j rjdtdx. 

We implicitly used the freedom of the form of the Lagrangian on the centre of to choose 
the kinetic term ^a? with a G M. Th en, because d = 0, we f ixed a = 1 and recovered the Euler- 


Poinc are equations (j2.9h . We refer to lMarsden et alJ 20071 1 i lGarcfa-Naranio and Vankerschaver 
2013l | for a similar construction. □ 


Provided the Legendre transformation is well-defined, the Euler-Poincare equation (|2.9p 
is equivalent to the Lie-Poisson equation (j2.6p . This can easily be seen by using the relation 
between the conjugate momentum L and velocity M 


^ ■= -w ^» 


or M := G g. 

oL 


2.2 Loop group and multi-time theory 


The idea of multi-ti mes for integrable system s was firs t intro duced by iDate et al.l 1983l | and 


further developed in iFlaschka et al 


Newelll 1985]. We will review here the key 


ingredients of this theory and then explain the links with the previous reduction theory. 


2.2.1 Loop groups and loop algebras 

The phase space is c onstructed from a particu lar inhnite dimensional Lie group, the polyno¬ 
mial loop group; see Pressley and Segal |l986| | for more details. For a semi-simple Lie group 
G, the associated loop group is G := Map(S'^,G), maps from the circle 5*^ with parameter A 
to group G. We will consider the elements of G through their Fourier series around A = 0, 
namely they will be polynomials with possibly an inhnite number of negative powers of A. 
The Lie algebra of G is then straightforward to construct. From the Lie algebra g of G the 
Lie algebra of G is g = Map(5'^,g). With the Killing form of the semi-simple Lie algebra g 
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one can also construct a pairing on g using the residue theorem. For two generic elements 

following calculation gives a simple form for the pairing 

I^{^^:VJ)X^^^dX = ( 2 - 11 ) 

i,j i,j i,j 

Notice that, the loop parameter A is different from the space variable x of the previous section. 
In fact, the finite dimensional Lie group G of the previous section will be replaced by the loop 
group G. We will thus use functions of A and x taking values in a finite dimensional Lie group 
G or Lie algebra 0 . 


2.2.2 Multi-time phase space 


Let the space-time manifold be a flat Riemannian manifold of countably infinite dimensions 
defined as N := lim„_^ool^"' and endowed with the standard metric gij = 5ij. The coordinates 
are denoted by a = (ai, 02 ,...). The choice of the letter a is to emphasize that there is no 
particular time or space variable, just coordinates on a manifold. In the following we will 
consider hyperplanes of N spanned by two variables, indexed by i and j. They are slices of 
the infinite dimensional space-time and they will be denoted Nij = span(ai,aj) C N. As we 
will see later, i is fixed at the beginning of the theory and j will be selected almost at the end 
to obtain a 1-1-1 PDF. In these slices there will thus be two interpretations for the physical 
meaning of a* and aj. We can either choose to set the spatial variable to be x := a,; and 
the time will be t := a-i or the reverse. The first interpretation, introduced by iDate et al. 


19^ : iFlaschka, et all |l983al lij: iNewelll |l98,^ . is common in the literature and leads to the 


standard Hamiltonian formalism. Surprisingly, the second interpretation seems to have never 
been notice d before in this p a rticular context, although it is an old idea in field theory; see 
for example ICaudrelierl |2015l ] ; ICaudrelier and Kundul |2015l | for other recent uses of this idea 
for the Sine-Gordon equation or for the inclusion of defects in integrable systems. We want 
to emphasize that the only difference lies in the fact that one of the coordinates is selected at 
a different stage in the theory. 

Going back to the full space-time manifold N, the phase space can in fact be understood 
in the context of c l assical field theory. Indeed, it is the first jet bundle J^{N,G). We refer to 


Gotav et all |l997t | ; iGastrillon Lopez et al.l |200ll ] for the general constructions of this space 
in field theory. In our case, the bundle structure N x G ^ N is trivial, thus the jet bundle is 
isomorphic to T*N®TG, the space of linear maps from TN to TG. With the left trivialization 
of TG, the reduced phase space is then {T*N ® TG)IG = T*N ( 8 > J. This phase space can be 
reduced by selecting a particular slice and is thus the phase space of a 1-1-1 PDF. A section 
of the bundle T*N (^q ^ N, namely a map M : A^ ^ T*N ( 8 ) 0 , corresponds to the projection 
of a map V : N T*N ^ TG only if the curvature of M vanishes for all a. The curvature 
is defined as the covariant exterior derivative of M with respect to M, viewed as a connection 
on this bundle structure and is given by = dM -|- [M, M]. In components, the curvature 
reads, for an element M = 


d^M = 


CXD 

^ A daj, 

*, 1=1 


( 2 . 12 ) 
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where ^a^ denotes the partial derivative with respect to Oj. The relation given by = 0 is 
called the zero curvature relation (ZCR) and contains an infinite number of constraints (one 
for each pair {i,j)) for the helds in M. 


2.2.3 Complete integrability 


Recall that the section M contains an infinite number of terms, each of them associated to 
a space-time direction, namely is associated to the direction Oj for every i. Each 
belongs to g and also has infinitely many components. The section M has therefore too 
many degrees of freedom, an inhnite number of infinite dimensional fields. For the complete 
integrability to arise, the number of independent fields of this system must be drastically 
reduced. It is done with the help of a very simple construction. We hrst dehne a particular 
loop algebra element with an essential singularity at A = 0, or an inhnite number of negative 
powers of A, by 

CXD 

:= ^Mi(a)A"\ (2.13) 

i=l 


where the Mi are sections of the bundle N x q ^ N. We will then use a shift operator 
and projections on the loop algebra. The shift operator is the multiplication by a power 
of A and the projections are projections on the two subalgebras with positive or strictly 
negative powers of A, denoted by P±. This decomposition in t wo subalgebra is crucial i n this 
construction and is at the root of the R-matri x forma l ism o f ISemenov-Tvan-Shanskiil 1983l | 


or the Marsden-Weinstein reduction exposed in iNewelll 
the other elements of the connection M as 


19851]. With these tools we can dehne 


mW := P+(A*M('^)). (2.14) 

With this particular construction, we have a one to one correspondence between M and 
and therefore the system does only depend on one loop algebra element We can thus 

expect that the inhnite number of constraints from the ZCR would be enough to sufficiently 
reduce the number of independent helds. What we actually expect to obtain from this ZCR 
is to reduce the inhnite number of Mj to only a hnite number where the freedom will be to 
choose the number of these independent helds by selecting a particular Once a z is 

selected, will be a function of through the implicit relations given by the ZCR 

()2.12p on the slice Woo = limj_).oo Wj- In this case they read 

= 0. (2.15) 


The formula holds true in a more general setting, and here is the proposition. 

Proposition 3. Choose a vector w € T^N, and, because (|2.12l) holds on every slice indepen¬ 
dently, it holds when one of the variable is Ooo := limj_).ooai. This can be written as a limit, 
where d^M : TN x TN ^ IS., 

lim =0. (2.16) 

n->-oo Y OCLn) 

The solution of this equation uniquely determines as a function ofM-w, the contraction 

of the one form M and the vector field w. In the case when w = ^ one obtains (|2.15l) . 
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Proof. By rewriting the ZCR using the definition of in the case of w = we have 

+ [P+(A*m(°°)),P+(A^m(°°))] = 0. (2.17) 

oaj oui 

Then, noticing that 

lim A"^P+(A^M(°°)) = 

J-fOO 


multiplying (j2.17l) by A ®, taking the limit j —)■ oo and together with A ■^P+(A®M(°°)) —)■ 0, 
we obtain the (j2.15p . □ 

The explicit computation of the Mj can be difficult depending on the Lie algebra g, i and 
the number of Mj that one wants to obtain. The Mi will also depend on through the 
derivatives, thus we leave the strict first jet bundle construction of this theory. One can 
also think of selecting a vector field w which could contains more elements and thus expect 
to obtain a higher dimensional integrable hierarchy. This is still an open problem because 
slices must first be extended to volumes and nothing is clear anymore. Up to this point we 
did not talk about dynamics of any of these fields. The ZCR (j2.15p will actually be the 
momentum-velocity relation needed for any dynamical interpretation. This will be done in 
the next section, depending on which formalism one wants to use. 

2.2.4 Prom ZCR to Lie-Poisson or Euler-Poincare equations 

In this section we will show that the Euler-Poincare or Lie-Poisson equations are the same as 
the ZCR after some preparatory steps. The first step is to select an integer k and compute the 
functions M^{M^^'>), Vj using the ZCR (j2.15p on the slice Nkoo- The second step is to select 
a n, thus to fix a slice where the 1-1-1 PDE will live. There is a third step before making 
any dynamical interpretation, namely decide what will be the space and time variables. There 
are only two choices and they will lead to two different formalisms: Hamiltonian if x := 
or Lagrangian if t := a^. 

Theorem 4. Let the connection M satisfy the ZCR = 0. For each slice Nij the restricted 
ZCR = 0, where + M^^'^daj, has an equivalent Hamiltonian or 

Lagrangian formulation: 

1. Hamiltonian: If the space variable is x := a,, for L := the ZCR (j2.15jl on 

the slice Nioo implicitly defines every Mj as a function of L. If the time variable is 
then t := aj, the Hamiltonian function is defined by hj{L) = f {L,M^^\L))dx and the 
associated Lie-Poisson equation is the ZCR on the slice Nij; see section l2.1.fl 

2. Lagrangian: If the time variable is t := aj, for M := M^^\ the ZCR (I2.15p on 
the slice Nj^o implicitly defines every Mi as a function of M. If the space variable is 
then X := a,, the Lagrangian function is defined by IfiM) = f (M, M^^^ (M))dx and the 
associated Euler-Poincare equation is the ZCR on the slice Nij; see section [2.1.3[ 
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This theorem relat es the well-known Hamiltonian formul a tion o f integrable hierarchies 


on loop algebras (see Newell 1985]; Semenov-Tvan-Shanskii 1983l |l to a new Lagrangian 


formulation through the generalized ZCR structure. Indeed, by using an abstract space-time 
manifold without any particular time or space variable, we were able to find the Lagrangian 
interpretation by avoiding an inverse Legendre transformation, which is the main obstacle to 
a comprehensive Lagrangian formalism of integrable systems. We will illustrate this theory 
in the next section but before that, we want to deeper address the question of the Legendre 
transformation. 


2.2.5 Legendre transformation 

We will briefly explain here how to relate the Lagrangian and Hamiltonian formalisms in our 
context with the Legendre transformation. In the Hamiltonian formalism the space variable 
X := Ofc is fixed once and for all and the hierarchy is then spanned with the time variable 
t := ttn- The conjugate velocity to the momentum L := is found by solving the recursive 
equations given by the ZCR (I2.15|] in the slice Nnoo-, up to the order n. For the standard 
hierarchies, n is always larger than k and the conjugate velocity will therefore contain up to 
n spatial derivatives. In the Lagrangian formalism, the time variable is fixed and the space 
variable has to be selected using the same procedure. Therefore, instead of spanning the 
integrable hierarchy with the time variable, it is spanned with the space variable. The crucial 
point is that in the Lagrangian setting, the equation will be written with more independent 
fields than it would have been in the Hamiltonian setting. Hopefully, the Euler-Poincare 
equations can be simplified into a single one for the momentum (or i equations if there are i 
momenta). This procedure of simplification can be seen as a Legendre transformation because 
the resulting simplified equation will be the same as if we started in the Hamiltonian side. 
Even if from a computational point of view the Lagrangian interpretation does not really 
differ from the standard approach, from a formal point of view it is of course very different 
and will be crucial in the development of the deformation of integrable hierarchies later in 
section [3l 


2.3 Application for standard hierarchies 

In this section we will be show how to use this formalism to recover integrable hierarchies 
such as the AKNS hierarchy. The main difficulty is to solve the ZCR or equivalently solving 
an implicit recursive system of equations. 


2.3.1 The AKNS hierarchy 


The Lie algebra for the AKNS hierarchy is s[(2) and the slices are given by x := oi for any 
other times t := Oj. For instance, selecting t := 02 will give the NLS flow and t := a flow 
which can be reduced to KdV/mKdV equation. We first recall the fundamentals of this Lie 
algebra. The basis matrices are 


0 1 
0 0 ’ 


/ 


0 0 
1 0 ’ 



0 

-1 ’ 


e = 


(2.18) 
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and the commutations relations 

[ej] = h, [e,h] = 2e, [f,h] = -2f. (2.19) 

We will also use the notation for the component of ^ in the Cartan subalgebra (which is 
an element A proportional to h) and ^ for the complement of 

Hamiltonian derivation In the Hamiltonian formalism one has to fix the space variable 
X = ai, thus the L operator is 

L = XMo + Ml. (2.20) 

By using the ZCR in the A^ioo slice we can express every other Mj, especially M := = 

A^Mq + A^Mi + AM 2 + M 3 in terms of Mi only because Mq will be constant. This amounts 
to solve the recursive relations defined by the ZCR (|2.15l) 

A° : d^Mo = 0 
A-^:a^Mi + [M 2 ,Mo] = 0 

A ^ : 9a:M 2 + [M 2 , Ml] + [M 3 , Mo] = 0 (2.21) 


A-' : d,Mi + [Mi, Ml] + [M,+i, Mo] = 0 


From the first equation Mq is a constant in ai and there are two choices for the value of 
Mq, an element proportional to the Cartan subalgebra or another element. The first, which 
is the most common in the literature, is related to the first grading of th e underlying K ac- 


Moody algebra, the second choice corresponds to the second grading; see Newell 19851 ] for 


more details linked to the Kac-Moody algebras. The two formulations are equivalent so we 
will stick to the first grading in this work and denote Mq = A where A = ih is then particular 
element proportional to the Cartan subalgebra. At this stage, one can remark that even if M 
is considered to be the independent variable, the ZCR (|2.15p still imposes some constraints on 
it. In this case the highest power must be constant. The third equation dxMi + [A, M 2 ] = 0 
also implies that Mi has no component in the Cartan subalgebra. M^ can then be computed 
using the fact that [A, [A,^]] = —for arbitrary ^ and reads M^ = —j[A,dxM^]. The 
parallel part of M 2 is found from the next equation by projecting out the perpendicular part 


dxM^ = -[Mi 


'Mi] = ^[[A,dxMi],Mi] = -^dx[Mi,[A,Mi]], 


where we used the Jacobi identity for the last step. With the vanishing (or periodic) boundary 
conditions, we have m| = — |[Mf*-, [A, Mf*-]]. From the very same equation. Mi can be 
calculated and is given by 

Mi = -^[A,9a;M2-^] - ^[A, [m|,Mi^]]. 

By denoting U := Mi, the ZCR on the slice A^i 3 , or equivalently the Lie-Poisson equation 
reads 


1 


1 


1 , 


dtU^ + - -dx[A, [U^, [U, [A, [/]]]] + -\U^, [U^, dxUi] = 0 


32 


( 2 . 22 ) 
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where the last two terms are in fact the same. This is a dynamical equation for the field 
U only and, by restricting its form, the KdV, mKdV or coupled KdV (cKdV) equations are 
recovered. Here is a summary, after rescaling time as t —>■ 4t, 


UKdV = 
UmKdV = 
UcKdV = 


0 U 
1 0 

0 u 
au 0 

0 u 
V 0 


lit QuUx T Uxxx — 0) (2.23) 

Ut Gcu Ux T '^xxx — 0? (2.24) 

Ut + 6uvUx + Uxxx = 0, Vt + 6uvVx + Vxxx = 0, (2.25) 


where cj = ±1 will give the focusing or defocussing mKdV. One can check that the hrst flow 
on Ni 2 is indeed the NLS equation (focusing or defocussing for a = ±1) is found with 


Unls 


0 u 
au 0 


iut + Uxx + au\ur = 0 . 


(2.26) 


The main reason why we only look at the 02 flow with the NLS reduction is that the other 
two reductions for mKdV/KdV are trivial. This can be seen in the calculation of the flow 
from the fact that one needs complex fields for a non-vanishing second flow 02 . This also 
explain the use of a complex constant A = i/i, in order to also include the NLS equation in 
the AKNS hierarchy. 


Lagrangian derivation In the Lagrangian formalism one has to hx t = 03 and then the 
M operator is 

M = \^A + + AH + IL (2.27) 

where A is still constant and U-^,V,W are three independent fields with values in the Lie 
algebra g and not the loop algebra. Solving the ZCR in order to find L = Mi is trivial and 
gives 

L = XA + U^. (2.28) 

Then, the Euler-Poincare equation, or ZCR in N 13 expands in three equations for the T part 

dtU^ - dxW^ + [U^, kPlI] = 0 

-dxV^+ [U^,V\\] + [A,W^]=0 (2.29) 

-dxU^ + [A, = 0, 

and two for the || part 

-dxW^^ + [U^,W^] = 0 

II , I (2-30) 

-dxV\\ + [U^,V^]=Q. 

One can easily check that this set of equations is equivalent to the KdV flow derived in 
the Hamiltonian formalism by expressing V, W in term of U only. This computation is the 
Legendre transformation from the velocity (C/, V, W) to the momentum [/, as described in 
section 12.2.51 
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2.3.2 50(3)-hierarchy 

This hierarchy, based on 50( 3), is not very w ell known in the literature and has recently 
been studied, for instance, bv iMal 20131 . 120141 ] . The Cartan subalgebra can be taken to be 
any element of the basis of the Lie algebra. We will choose 63 , where the basis of so (3) is the 
standard ( 61 , 62 , 63 ) with the commutation relations 


[ 61,621 = 63 , [ 61 , 63 ] = 62 and [ 62 , 63 ] = 61 . 


(2.31) 


Following the AKNS scheme, we use A = ies, where 63 is taken as the Cartan subalgebra 
basis vector. 


First flow of the hierarchy The Euler-Poincare equation is found using the usual two 
elements 


and reads 


L = AA + U^, M = A^A + AU^ + V, 

-d,u^ + lA,v^] = o 
-d^v'' + iu^,v^] = o. 

After computing the Legendre transformation, V can be expressed as 


v^ = -[d,u^,A], 


V\^ = -^[U^,[U^,A]]. 


The 1+1 PDE, when U = uei + U 62 , finally reads 


1 


idtu = Uxx + ^(«^ + u^)u. 


(2.32) 


(2.33) 


(2.34) 


(2.35) 


This equation is a modification of the nonlinear Schrodinger equation where jtt]^ is replaced by 
the difference Re(u)^ — Im(n)^ together with appropriate conjugations. This equation seems 
to be new, but does not have the f7(l) phase symmetry. The other cho ice of A, riamely A = 63 
and U = uei + U 62 leads to two coupled equations already derived in Ma 20131 . 2014 j. 


Second flow of the hierarchy The second flow has the now usual elements 

L = AA + U^, and M = A^A +A^U^ + AV + W. 

The Euler-Poincare equation is then 

dtU^ - dxW^ + [U^, kpll] = 0, -dxWW + [C/^, W^] = 0 

-dxV^ + [t/^, pii] + [A, W^] = 0, -dxVW + [U^, E^] = 0 

-dxU^ + [A, E^] = 0. 


(2.36) 


(2.37) 
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After the Legendre transformation, we obtain 

A], yii = [u^,A]] 

= dlu^ + ^[A, [U^, [U^, [U^,A]]]], LL" = [U^,d,u^]. 
With U = uei + ve 2 , the coupled equations for u and v read 


(2.38) 


dt^U = Uxxx — 2 ^^^^ — — 2uVVx 

dt^V — Vxxx A ^^x ‘2uVUx' 


(2.39) 


This equation is then a coupled mKdV equation for u and u. also found in iMal 201311 an d 


previously in the study of coupled modified KdV equation of iTsuchida and Wadati 
One can further simplify this equation by setting u = v and recover the modified KdV 
equation. This illustrates the fact that the first members of different hierarchies are sometimes 
the same, owing in this case to the isomorphism between the Lie algebras s[(2) and so(3). 


3 Deformations of integrable hierarchies 


Apart from classical completely integrable syste ms, there exist other intere sting systems such 
as the Camassa-Hohn equation, first derived in Camassa and Holm 1993| | and which admits 
peaked solitons, or peakons. It is of common agreement that their complete integrability is of 
different flavour than classical integrable systems. From our point of view, the main difference 
is that they d o not f i t into the present loop group approach. Notice that a successful try had 
been made bv ISchifll 19981 ] where he mapped the CH equation to a negative flow of the KdV 
hierarchy using a reciprocal transformation. Indeed, the ZCR of the CH equation does not 
have a constant element for the highest power of A , but the dynamical field itself as the ZCR 
is given by (see for example Hone and Wang 20^) 


L = 


-m\ + j 0 


M = 


2^x 

umX + ju— 


-u-\X 1 


(3.1) 


Notice that we used the momentum m = u — a^Uxx where is the length scale parameter. 
By letting —>■ 0, the dispersive CH equation reduces to the KdV equation but the ZCR of 
the dispersive CH equation will not converge to the ZCR of the KdV equation, written in the 
AKNS matrix form. This means that the CH equation has a different integrability flavour 
than the KdV equation. On top of that, their is no AKNS hierarchy for CH type equations 
despite their close relationship with the AKNS hierarchy. These differences are important 
because, from the equation standpoint, CH is a deformation of KdV, but from the integrable 
system theory, they seem to have nothing to do with each other. We will hereafter show how 
to deform the integrable system theory developed above such that the CH equation, amongst 
others, can be recovered. This will lead us to a definition of a weak integrability in the next 
section. 
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3.1 Weak integrability 


The present formulation of integrable systems using Lagrangian mechanics is the key ingre¬ 
dient for a theory of deformation of integrable systems. By deformation, we mean replacing 
the LP' norm by the norm. This procedure introduces a length scale parameter a such 
that when a —)■ 0 the norm becomes the Lp norm. We first recall that, in classical me- 
ch anics, the metric can be d efined by the Lagrangian if it is in a quadratic form. We refer 

proper 


to iMarsden and RatiJ [l999^ for a precise account of this fact. We do not have a 


quadratic form here because the momentum is related to the velocity through the compli¬ 
cated recursive relations (I2.15P as shown before. Nevertheless it is the best place to introduce 
the norm. The deformed Lagrangian can therefore be defined with the Sobolev norm as 


= j{M^^\ (1 - a^dl)M^\M^^^))dx. 


(3.2) 


For the last equality we rewrote the norm using an integration by parts and vanishing bound¬ 
ary conditions. This exhibits the Helmholtz operator A := 1—that we will use throughout 
the rest of this work. This type of deformation is common in the literature and allows singular 
solutions to exist. Indeed, the Green’s function of A, given by l{2a) is is the f amou s 


peakon solution of the di spersionless CH equation. We refer to Camassa and Holm 1993| | 


Holm and Marsden 20051 ] and subsequent works for different aspects of peakon solutions. 


After having modified the Lagrangian with the Sobolev norm, the standard Euler-Poincare 
equation (|2.9I) can be derived and reads 


Vi, 


(3.3) 


for A := 1 — oPd"^. We want to emphasize that M^) had already been fixed and that the ZCR 
(|2.15[) on A(joo had been used to express the corresponding momentum before 

the Euler-Poincare equation was derived. This modification does not change the implicit 
relation but changes only the Euler-Poincare equation. From this modified Euler- 

Poincare equation (|3.3p . one can try to compute the Legendre transformation in order to, 
for instance, recover the CH equation. The deformations will be computed in section 13.21 
but first, one has to be careful with the integrability property or, saying differently how to 
come back to the ZCR interpretation from the Euler-Poincare equations (|3.3p . Indeed, this 
is not a trivial operation and it can only be done after defining a weak ZCR, or projected 
ZCR. For convenience here, we will only work with slices of the form Nn where x := ai as 
fiows with higher order space variables are more complicated and not much used for physical 
applications, except maybe for the derivative NLS equation with t = 04 and x = 02 ; see 
Flaschka et al. 1983a j: Kaup and Newelj 1978| |. 


Definition 5. LetVfj. : g —)• g 6e a projection operator for polynomial loop algebras depending 
only on A, the element in the Cartan subalgebra of g and k, which defines a slice Ni^ for the 
equation in considerations. For an arbitrary Z = X)i=-cx) ^ 0? projection is given by 


Pffc(Z) = Z - - \^-^P+{\^-'^{Z,A)A/{A,A)) 


k-i- 


a-fc/ 


(3.4) 


where P+ stands for the projection onto positive powers of X. 









































17 


This projection corresponds to removing the Cartan subalgebra element of the Z^-i and 
the full element Z^- We can then naturally define the notion of weak integrability. 

Definition 6. The Lie algebra value two form + [M^^\ M^^'>])ds/\dt is said 

to be a weak ZCR if its projection under vanishes, namely 

= 0, (3.5) 


or equivalently 


= 0. (3.6) 

If the weak ZCR is equivalent to a PDE, the PDE is said to be weakly integrable. 

From this definition, we see that apart from being proportional to a power of a the 
projected terms are proportional to the two highest powers in A that appear in the ZCR. 
The equations proportional to these powers of A are always naturally verified if a = 0, 
indeed, we will see that one of them is of the form [U, U] = 0, for a matrix U. Notice 
that the weak ZCR (j3.5p is still written in term of the velocities in the Lagrangian, but 
can be Legendre transformed to be expressed with the momentum fields only. After the 
Legendre transformation, the ZCR will still be a weak ZCR but it will be equivalent to some 
nonlinear PDE, as the CH equation. In order to completely understand the integrability of 
these systems, the isospectral problem must be understood. In the non-deformed case, the 
isospectral problem is standard but after the deformation its correct formulation is still an 
open problem. 

We will display the exact projected terms for each equation later, but first we want to 
sketch a possible research direction for this problem. For simplicity, we denote by L and M 
the two loop algebra elements of the ZCR. The trick is similar for the Manakov triple used 
for 2 + 1 integrable PDEs, namely instead of projecting out some terms of the ZCR, one can 
recast it inside the commutator as 

ALt -Mx + [AT, M - X^L] = 0, (3.7) 

if k corresponds to the flow of M = Then, a direct calculation gives the spectral 

problem for a wavefunction t, A) 


V’t = MV’ — X^Lip 

= alV'- 

This spectral problem is then isospectral if and only if 

{L4>)x = {Lx + LAL)V’ = 0. 


(3.8) 


It remains an open problem to interpret and solve this spectral problem as it is rather different 
from the Manakov triple. We will not address this problem here but rather focus on the 
geometrical interpretation of the deformations of integrable systems. 
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3.2 Deformed AKNS hierarchy 

Following the derivation of the AKNS hierarchy done in section r2.3.1l but with the deformed 
Lagrangian we will derive the Camassa-Holm equation as well as other equations. In the 
Lagrangian formalism one has to fix f := 03 and the M operator is then 

M = \^A + + AK + VF, (3.9) 

where A is still constant and !/■*■, V, W are three independent fields. Solving the ZCR in order 
to find L = is trivial and, after applying the A := 1 — operator to L, we readily 
have 


L := AL = XA + AU-^. (3.10) 

The Euler-Poincare equation, or ZCR in the slice iVi 3 , expands in four equations for the T 
part 

A° : dtAU^ - + [AU^, ipll] = 0 

Ai : + [AK^, Pll] + [A, IF^] = 0 ^ 

+ [A,V^] = 0 
X^ : [A,AU^] + [U^,A]=0 

and three for the || part 

A° ; -5,,if" + [AU^,W^] = 0 

A^ : + [AU^, = 0 (3.12) 

A^ : [AU^,U^] = 0 

where the last equation of both systems are no more trivially satisfied and has to be projected 
out with the projection operator Indeed, one can check that the projection exactly 

removes these two terms. This example illustrates the fact that the lack of complete integra- 
bility is rather small for the deformed equations and that with = 0, the projection does 
nothing. 

The Legendre transformation can now be computed to obtain equation such as the CH, 
mCH and the new CH-NLS equation. First, the A^ equation of (13.lip gives 

= -^[AdxU^], 

and then the A equation of (j3.12p yields 

dxVW = [AU^,V^] = -^[AU^,[A,dxU^]]. 

This equation can only be weakly solved as 

v\\ = -^d-\[AU^,[A,dxU^]]). 










19 


W~^ is non-local 

= -\[A,a^V^\ - 1[A, [ylUf/-L|] = -\8lu^ + i[A, [8j‘(|Af/^,|A,a,(7-L|]),Af/|], 
as well as the parallel part of W, which reads 

d,w\\ = [W^,W^] = -l[W^,dlU^] - [A, [d-\[KU^, [A,d^U^]]),AU]]]. 

3.2.1 CH and mCH equations 

In order to obtain the standard form of integrable wave equations we have to fix the form 
of U and rescale the time t At. For UxdV defined in (j2.23jl one can check that the CH 
equation is recovered 


mt + 2mux + ruxU + Uxxx = 0 , 


(3.13) 


with the Uxxx dispersive term; see Camassa and Holm 1993l |: Dullin et al. 2004]. The corre¬ 
sponding weak ZCR after simplification of the complex numbers is given by 


L = 


A m 
1 -A 


M = 


A^ - \Xu - \ux \^\ux - \mu + \uxx 
A2 - in 


—A^ -|- iAn -|- 4 


4 


(3.14) 


Notice that the term that we need to project out in order to recover the CH equation (|3.13p 
from the previous weak ZCR is 


-A2 2A3 
0 A^ 


A2([H,AH] + [U,A]) + X^[AU,U] = a\xx 

For UmKdV: defined in (12.241) . the dispersive mCH equation is obtained 

mt 2cT[m(n^ - a^ul)]x + Uxxx = 0. 


(3.15) 


In this case, projects only [AC/-*-, C/-*-] = 0, thus the weak ZCR has only a projection for 
the A^ term. The M operator reads 


M = 


A^ — ^A(n^ — a^u^) 

X^u — \Xux — \m{u^ — o?u^) + 4 


X^u + X\ux — \m{u^ — cP'u^) -|- 4n; 


Ux 


—A^ -|- 4 A(u^ — oP'u^x) 


4 ^xx 


Similarly to the CH equation, the term that we need to project out in order to recover the 
dispersive mCH equation (I3.15P from the previous weak ZCR is 


A2([A, AU] + [C/, A]) + X^[AU, U] = a^ux 


-X^ 2X^' 

-2A3 A2 


The mCH equation is already known t o be i ntegrable with a linear disp ersion see iQiao 


20071 ] for a recent derivation. We refer to lOiaol [20071 ]: iFokas and Liul |199 (tI ] and references 


therein for more details on this equation, we will not investigate it further in the present work. 
The mCH equation with a third order dispersion term as we derived can easily be related to 
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the linear dispersion by a change of variable x ^ x + ct for an appropriate value of c. Finally, 
the general UcKcLV defined in (I2.25P gives the coupled mCH equations 


(3.16) 


rtit + 2[m{uv — a^UxVx)\x — 2m{uvx — Uxv) + Uxxx = 0 
nt + 2[n{uv - a^UxVx)\x + 2n{uvx - Uxv) + Vxxx = 0. 

These coupled equations have recently been found and studied by Xia and Qiaol 2015]. The 
weak ZCR can be calculated but we will not display it here. 


3.2.2 Deformation of the NLS equation 


The deformation of the NLS equation, the first flow in the AKNS hierarchy, can be computed 
and will give a new weakly integrable equation that we will call the CH-NLS equation. Using 
the previous calculations of M = \^A + At/-*- + V and the NLS form of U 


^Unls 


0 m 
am 0 ’ 


(3.17) 


for complex valued u and a = ±1 for the focusing or defocussing case, we obtain the CH-NLS 
equation on the slice N 12 

imt + Uxx + 2am{\u\‘^ — a‘^\ux\‘^) = 0, m = u — a^Uxx, a = ±1. (3.18) 


In term of u only it is given by 


9 I 19 9 I 

iut — ia Uxxt + Uxx + 2au\u\ — 2aa u|tta 


— 2aa^Uxx\u\^ + 2aa^Uxx\u3 


= 0 


13.191 


but is not an evolutionary equation for u, as all the other deformed equations. The weak ZCR 
of the CH-NLS equation is 



m 


m 

-iX ’ 


M = 


i\^ + §(|u|^ — a^\ux\^) 


Xu + ^Ux 


Xu — ^Ux 
—iX^ — §(|up — 



(3.20) 


and, in contrary to the mCH equation, the projection with respect to the term remains, 
because of the complex valued fields. This equation is also Hamiltonian with its Hamiltonian 
structure given by the non-canonical NLS Hamiltonian structure 


J = 


2amd^ 

dx + 2amd~^m 


dx + 2amd^ ^ m 
2amdx^m 


and its associated Hamiltonian 


P = i 



mux)dx. 


(3.21) 


(3.22) 


The Hamiltonian has an interpretation of momentum for the field m, standard in the theory 
of the NLS equation. This Hamiltonian structure is actually the same as for the NLS equation 
and the modified definition of the momentum P leads to the CH-NLS equation instead of the 
NLS equation. Note that the mass M = J \m\‘^dx is also conserved and could be associated 
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to the symmetry of the CH-NLS equation. The interpretation of a mass and momentum 
for M and P is not clear, as the Hamiltonian structure J does not produce the flow of space 
translations mt = mx and phase shifts mt = im. Notice that there is no Galilean symmetry 
for this equation. This is linked to the non-integrability of the CH-NLS equation. Indeed, 
despite the weak integrability, this equation seems not to be completely integrable, as a second 
compatible Hamiltonian structure as well as its associated Hamiltonian are missing. For the 
NTS equation, the second Hamiltonian structure, which is a canonical Hamiltonian structure, 
would generate the symmetry associated to the ma ss M and mome ntum P. We will not 
investigate this equation further here, but we refer to Arnaudon 2016| | for more details about 
this equation. 


3.3 Deformation of SO{3) hierarchy 

Following the same procedure as for the deformation of the AKNS hierarchy, we proceed with 
the 50(3)-hierarchy. 


First flow of the hierarchy The L and M element are 

L = XA + AU^, M = X^A + XU^ + V (3.23) 

and the corresponding Euler-Poincare equation is 

dtAU^-dxV^ + [AU^,v\\] = 0 
-dxU^ + [A, = 0, + [AU^, P^] = 0. 

Then the Legendre transformation gives 

P^ = -[dxU^, A], pll = -d-^[AU^, [dxU^, A]]. 

The deformation of S03-NLS equation (I2.35P then reads 

• _— I f / 2 2 2 I —2 2 —2 \ / o o \ 

imt = Uxx + [u — a Ux + u —a Ux) ■ ( 3 . 25 ) 

Even if this equation seems to be new, the lack of 17(1) symmetry makes it less physically 
relevant we thus leave the analysis of this equation for future works. 


Second flow of the hierarchy For this flow, the L and M elements are 

L = AA + AU^, L = X^A + X'^U^ + AP + IP (3.26) 

and the associated Euler-Poincare equation is 

dtAU^ - dxW^ + [P^IPII] = 0 , -dxW\\ + [AU^W^] = 0 

-dxV^ + [At/^, pll] + [A, IP^j = 0, -d,,pll + [AP^P^j = 0 

-9,t/^ + [A,P^] =0. 


(3.27) 











22 


After computing the Legendre transformation we obtain 

= -[d,U^,A], = -d-^[AU^, [d,U^,A]] 

= dlU^ + [A, [AU^, [d-\AU^, [d,U^,A])]]], fL" = d-^[AU^,dlU^]. 

The equation for AU = mei + ne 2 is then 

mt + [miu^m + v^n)]^ - {uv^ - vux)v + u^xx = 0 
nt + [n{uxm + Vxn)\x + {uvx - vux)u + Vxxx = 0 . 

This equation is similar to (13.161) except for the third term; there might exist a transformation 
between the two. If one restricts the form of U by setting v = u, the equation becomes 

mt + Uxxx + 2m?Ux + mx{u^ — = 0 (3.29) 

which is the modified CH equation (I3.15P . This result is compatible with the classical hierarchy 
which gave the modified KdV equation. As in the classical case, the difference in terms of 
the form of the equations between the two hierarchies arises when the full UcKdV element is 
considered. 


3.4 Limiting case: —)■ oo 


The limit —>■ oo is also interesting and corresponds to the high frequency limit when 

e ^ oo in the change of variables x —>■ ex,t et. For the Camassa-Holm equation i t gives 
the Hunter-Saxton (HS) equation iHunter and SaxtonI |l99l[ |: iHunter and Zhend [1994 1 


Utx + UUxx + -Ux = 0 - 


(3.30) 


This limit corresponds to modifying the Sobolev norm with its equivalent norm f u^dx. In 
the case of cubic equation, the high frequency limit which corresponds to this pairing has a 
different form, namely x —>■ ex, t —>■ e^t. The limit of CH-NLS equation (|3.18l) reads, after 
rescaling the time t —> 2at 


^^xxt — ^xxl'^x 


2 


(3.31) 


This equation seems to be new, but its analysis is beyond the scope of this work. The same 
limit in the mCH equation yields 


'^xxt — '^xx'^x 


or 


'^xt — 3^^ 


or 


Vt = -V 


(3.32) 


which is now a ordinary differential equation for v = Ux, thus not interesting for us here. 
The high frequency limit for the other equations will not be displayed here, but might be 
interesting for further studies. 
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4 Conclusion 


In this paper, we deformed the classical integrable system theory in order to derive nonlinear 
equations with nonlocal terms such as the Camassa-Holm equation. The central point of this 
deformation is in a systematic insertion of the Sobolev norm in hierarchies of integrable 
systems. In order to achieve this goal in a systematic way, we replaced the standard 
norm by the Sobolev norm in a Lagrangian which should describe an integrable hierarchy. 
In order to derive such a Lagrangian formulation of integrable systems, we first came back 
to the roots of integrability, written in term of multi-times and loop groups. Then, we took 
a slightly different direction by forgetting for a moment the time or space interpretation of 
the coordinates in the so-called multi-time space, associated with the different flows of the 
hierarchy. This step allowed us to reverse the choice of time and space in the construction 
of the hierarchy of equations and thus to interpret the usual zero curvature relation as an 
Euler-Poincare equation. This Euler-Poincare equation is rather special as it is written on the 
dual of a loop algebra and with an extra term coming from a derivative cocycle, responsible 
for the spacial derivatives in the resulting nonlinear equations. This exchange of time and 
space in the derivation of the hierarchy produced equations with more dependent fields, and 
thus led to coupled equations rather than a single equation. These coupled equations can be 
simplified by expressing the extra fields in term of a single one (in the case x = oi) and is 
understood as a Legendre transformation. This Legendre transformation then recovers the 
standard ZCR with its Hamiltonian Lie-Poisson interpretation. Notice that we treated here 
the simplest case of quadratic Lagrangians, written on se mi-simple Lie algebras. Extensions 
such as non semisimple Lie algebras treated for example in iMal 2009(1. or sem i -direc t products 
extensions, studied in the context of the CH2 equation by Holm and Ivanov 201Cl| |. could be 
interesting directions to explore in this context. 


Erom this viewpoint, the best choice that we have for the inclusion of the norm is 
in the Lagrangian. Then the Legendre transform will give the usual notion for of nonlinear 
equations with a single field, and a notion of deformation of hierarchies of integrable systems. 
The obtained equations do not form a hierarchy by themselves as some of them will not 
be integrable. The complete integrability in term of the standard ZCR is thus altered for 
this deformation of hierarchies. Indeed, the classical ZCR is no more valid but a notion of 
projected ZCR can be defined such that the equivalence between the deformed PDE and 
the deformation of the ZCR is retained. These projected terms are always proportional to 
the highest powers in the spectral parameter A, as well as to a^. Erom here, the link with 
an isospectral problem seems to be lost, but we suggested a possibility to incorporate this 
extra terms in the commutator of the ZCR, provided some extra conditions on the spectral 
problem hold. Even if most of the equations found in this work are already known, this 
deformation theory relates more closely the deformed equations with their classical limits 
given by —)■ 0. This systematic approach led us to a classification in term of deformed 
flows of deformed equations summarized in table dLH). In this classification, a new equation 
called the CH-NLS equation (II.ip was found as a deformation of the NLS equation. The 
CH-NLS equation could not yet be shown to be completely integrable and therefore raises the 
question of understanding the link between the present notion of weak integrability and the 
standard complete integrability. A final comment regarding the deformed equations is that 
there might be a possibility of using the KAM theory of PDEs for our deformed equations, 
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such that the one described in Kuksin 2000l |. The main problem here will be to treat the 


nonlocalities arising from the use of the inverse Helmholtz operator, as we should view m 
as the dynamical variable for applying this theory. A possibility to overcome this would 
be to expand the nonlocal terms up to some order in a and restricting the validity of the 
approximated equation for solutions with low wavenumbers compared to the scale given by a 
. Then, maybe the KAM theory could give some insights into the approximated equations, 
such as deriving approximated solutions for them. 
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